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Abstract. An n-axy fe-radius sequence is a finite sequence of elements taken 
from an alphabet of size n such that any two distinct elements of the alphabet 
occur within distance k of each other somewhere in the sequence. These se- 
quences were introduced by Jaromczyk and Lone to model a caching strategy 
for computing certain functions on large data sets such as medical images. Let 
/fe(n) be the shortest length of any fc-radius sequence. We improve on earlier 
estimates for /fc(n,) by using tilings and logarithms. The main result is that 
fkiji-) ~ s (2) as n — > 00 whenever there exists a tiling of by a certain 

cluster of k hypercubes. In particular this result holds for infinitely many k, 
including all k < 194 and all k such that fc + 1 or 2A; + 1 is prime. For certain 
k, in particular when 2A; + 1 is prime, we get a sharper error term using the 
theory of logarithms. 



1. Introduction 

Let k and n be positive integers. An n-ary k-radius sequence is a finite sequence 
ao, oi, . . . , Um-i of elements taken from an alphabet F = {0, 1, . . . , n — 1} with the 
following property. For all distinct x,y G F, there exist i,j E {0, 1, . . . , m — 1} such 
that tti — X, Uj ~ y and |i — j| < k. In other words, any two distinct elements of 
F occur within a distance of k of each other somewhere in the sequence. We may 
think of a window of length fc + 1 sliding along the sequence; the condition of being 
a fc-radius sequence says that we will always be able to find a window containing 
both X and y, whatever our choice of x and y may be. For example, a 5-ary 2-radius 
sequence of length 7 is: 

0, 1, 2, 3, 4, 0, 1. 

We write fk{n) for the shortest length m of an n-ary fc-radius sequence. The 
example above shows that /2(5) < 7 (and in fact /2(5) = 7). 

The notion of a fc-radius sequence was introduced by Jaromczyk and Lone 
Such sequences describe a First-In First-Out caching strategy for computing func- 
tions that require computations on all pairs taken from a set of n large objects 
(such as medical images), where at most k + 1 objects are cached at any one time. 
We believe that the problem of constructing short fc-radius sequences is interesting 
as pure combinatorics; this paper establishes surprising (to us) connections to the 
theory of tilings and logarithms, and to properties of cyclotomic fields. 

Let k be fixed. This paper is concerned with establishing the asymptotic growth 
of fk{n). It is obvious that /fc(n) < 2(2), as the concatenation of all length 2 
sequences x, y with 0<x<2/<nisa fc-radius sequence. Moreover (see Section [2| 
it is not difhcult to show that \{^^ < fk{n). So fk{n) seems to grow like (2), and 
we would like to determine the limit lim„_i.oo fk (n) / (2) (if it exists) . 
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Ghosh ^ proved that 

(1-1) /iH = < 

Jaromczyk and Lone [TT] showed that 



1 when n is odd, 



n/2 when n is even. 



Chee, Ling, Tan and Zhang [2J observed that the exphcit upper bound derived by 
Jaromczyk and Lone is poor for small parameter sets, and provide good construc- 
tions of 2-radius sequences when n is small. We modify the argument of Lone and 
Jaromczyk to improve the error term in their estimates (see Section [4| , showing 
that 



(1-2) 



/2(") 



+ 0(n^-^^^). 



Moreover, we show that there are infinitely many values of n such that /2(f^) 
iQ)+0(n). 

Turning to the case when k > 2, Jaromczyk and Lone [11] show that 



/fc(n) < 



1 



2[(fc + l)/3j V2 



(1+0(1)). 



We show that 
(1.3) 



whenever there exists a Z-tiling of Z"'*^' by copies of a certain fixed cluster of k 



hypercubes (Theorem 6.2). In particular, the result holds for all k < 204 except 
possibly k = 195, and the result holds whenever at least one of fc + 1 or 2fc + 1 is 
prime. For many of these values of k (including all values of k with 2fc+ 1 prime, or 
with fc + 1 prime and 8|fc) we are able to improve the error term in ( 1.3 ); moreover 



0(n). 



for these values of k we show that there are infinitely many values n such that 

(1.4) 

The structure of the remainder of the paper is as follows. In Sj2] we establish 
simple lower bounds for /fe(n), and recall some of the number-theoretic tools that 
we shall need later. We then describe a construction method (generalising that of 
[llj ) that for a prime number p produces p-ary fc-radius sequences via coverings 
of Z* by arithmetic progressions. We use this to produce a tight upper bound for 
f2{p) (Theorem [XT]), and deduce ( [Q] ). 

In Sjsjwe recall and develop the theory of logarithms (bijections from {1, . . . , fc} 
to Zfe such that f{ab) = /(a) -1-/(6) whenever all three function values are defined). 
We introduce the notion of a special KM-logarithm which imposes precisely the right 
parity conditions that we need for a later construction of fc-radius sequences. We 



explain the relationship between logarithms and certain tilings of Z*" (Theorem 5.3 ). 
We prove our main result (Theorem 6.2 1 in Sj6j that the existence of a suitable tiling 
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implies the estimate ( 1.3 1; in particular, the existence of a logarithm of length k is 
enough for this estimate to hold. 

In ^we define the notion of a /c-radius prime, and show that if n is a fc-radius 



prime then we get the estimate (1.4 1 for fk{n) with best-possible main term and a 
good error term (Proposition 7.1). We establish a link with logarithms by showing 
that fc-radius primes exist precisely when there is a special KM- logarithm of length fc 
(Theorem 7.2 1. In ^jsjwe derive the density of fc-radius primes for any fixed fc, and 



use this result to improve the error term in our estimate (1.3) for those values 
of fc where a special KM-logarithm of length fc exists. In f|9 we illustrate the 
applicability of our constructions by computing lists of all fc up to 300 for which 
there exist logarithms or special KM-logarithms. Finally, in ^10\ we discuss some 
of the various open problems that remain. 

2. Preliminaries 

This section contains some elementary lemmas on fc-radius sequences, and briefly 
describes some of the notation and results from number theory we will use. 

2.1. Elementary results on fc-radius sequences. 

Lemma 2.1. Let k, n and x be positive integers. Then 

(i) /fc+i(n) < /fc(n); 

(ii) fk{n) < fk{n + x). 

Indeed, the inequality (ii) may be improved to 

fk{n) < fk{n + x)~ \{x/{n + x))fk{n + x)] < {l ^ x/ {n + x)) fk{n + x). 

Proof. The first inequality of the lemma follows since every fc-radius sequence is a 
(fc -I- 1) -radius sequence. 

To establish the final inequality of the lemma, note that we may construct an 
n-ary fc-radius sequence by taking the shortest (n + x)-ary fc-radius sequence and 
deleting every occurrence of the x most common of the n + x alphabet symbols. 
The X most common symbols together occur at least \ {x / (n + x)) fk{n + x)~\ times, 
and so the final statement of the lemma follows. □ 



Lemma 2.2 (Jaromczyk and Lone, [11 ). For any positive integers fc and n, 

< fk{n). 



1 /n 
fcU 

Proof. Let ai, a2, . ■ ■ , am be a fc-radius sequence. The fc-radius condition implies 
that every unordered pair of distinct alphabet symbols must occur at positions i 
and j in the sequence, where I < i < j < m and \i — j\ < fc. The number of 
positions i,j where 1 < i < j < m and |j — j| < fc is strictly less than km, since 
there are m—1 choices for i and once i is fixed there at most fc choices for j . Hence 
km > (2) , as required. □ 

We remark that Jaromzcyk and Lone [llj show that the linear term of this bound 
can be improved, by considering the elements at the start and end of the sequence 
in more detail; they use this technique to give an explicit improved lower bound in 
the case when k — 2. 
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2.2. Some number theory. We use various pieces of standard terminology and 
notation from number theory. We use Z„ to denote the ring of integers modulo n, 
represented by the elements 0, 1, . . . , n — 1. The multiplicative group of invertible 
elements we denote Z* , its order being given by the Euler totient function (p{n). 

The letter p will always indicate a prime number. We use for the Legendre 
symbol. A natural number n is said to be y-smooth if each prime p dividing n 
satisfies p < y. The function 7r(a) counts the number of primes in the interval [1, a], 
and uj{n) denotes the number of distinct prime factors of n. 

The following theorem will be useful in moving from estimates valid for primes p 
to general estimates for natural numbers n. 

Theorem 2.3 (Harman, [9 , Theorem 10.8). Let a and s be coprime positive inte- 
gers. Then for all sufficiently large n, there exists a prime p such that p = a mod s 
and n < p < n + n"'^^^. 

3. A CONSTRUCTION 

In this section, we describe a construction method that we shall use to derive our 
asymptotic results. The method is a straightforward generalisation of a construction 
of Jaromczyk and Lone 111]. 

Let p be a prime number. For an element d e Z*, we define the set Bk p{d) C Z* 

by 

Bk,p{d)=d{±l,±2,...,±k) 

= {-kd,-{k - l)d, -2d, -rf, rf, 2d, . . . , (fc - l)d, kd}. 

For any d G Z*, let be the periodic sequence (0, d, 2d, . . . ,{p ~ l)d, 0,d, . . .) 
of period p. Note that distinct elements x,y £ Zj, appear at distance k or less 
somewhere within if and only if y — x € Bk^p{d). 

Theorem 3.1. Let _D C Z* have the property that the sets Bk^p{d) where d (z D 
cover Z* . Then there exists a p-ary k-radius sequence of length \D\{p + fc — 1) + 1. 

Proof. Let D = {di,d2, . . . ,du}, where u = \D\. Let be a finite sequence of 
length p -\- k consisting of p + k consecutive terms of the periodic sequence s^. . For 
i > 2, we choose the first element of to be equal to the final element of ti_i. Note 
that the sequence contains all pairs {x,y) with y ~ x £ Bk^p{di). The sequence 
tit2 • • • t„ of length u{p+k) produced by concatenating the sequences is therefore 
a fc-radius sequence, by the 'covering' property of D. But the final element of each 
sequence is equal to the first element of the sequence t^+i, by our choice of the 
shift of Sfj. for i > 2. We may always replace any repeats in a fc-radius sequence by 
single occurrences, to produce a shorter fc-radius sequence. There are u — 1 places 
where we may do this in our construction, so we have a fc-radius sequence of length 
u{p + fc — 1) + 1, as required. □ 

4. 2-RADIUS SEQUENCES 

In this section, we establish bounds on f2{n). Our methods are similar to those 
of Jaromczyk and Lone [11) , but our slightly different approach gives a better error 
term than theirs. This section is a good introduction to our more general methods 
in later sections of the paper. 
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Theorem 4.1. Let p be an odd prime number, and let £ be the multiplicative order 
of 2 modulo p. Then 



f2{p) < { 



+ + l tf its odd; 
i(p+l)(p-l)(l + f) + l 1/^ = 2 mod 4; 
+ 1) + 1 i/^ EE mod 4. 



Proof. Suppose that £ is odd. Then the theorem follows by Theorem |3.1[ provided 
that we can find a set of {p — 1)(1 + l/£)/4 sets B2.p{d) that cover Z*. Let (2) be 
the multiplicative group generated by 2 in Z*. Now Z* is the disjoint union of the 
cosets of (2): 

(4.1) Z; = CiUC2U---UCt 

where Ci = {ci,2ci, . . . ,2^~^Ci} for some elements q € Z*. Since |Ci| = £ for 
aU z e {1,2, . . . we find that t = (p - \)/£. For any i £ {1,2, . . . ,<}, we 
have that —d = Cj for some j € {1,2, .. . ,t}. Since £ is odd, —1 ^ (2) and so 
i j. So, by relabelling if necessary, we may assume that C2S+1 = ^C2s for all 
sg {l,2,...,t/2}. Define 

D ^ Di\JD2\J---\JDt,2, 

where 

D,^{c2s2'''\Q<i<{£+l)/2) 

for all s G {1,2,..., t/2}. Since d, 2d £ i?2,p(d) for any d e Z*, we see that the 
blocks i?2,p(rf) with d £ Dg cover C2S. Since —d,—2d G B2^p{d), we see that the 
sets B2.p{d) with d £ Dg cover — C2S = (^28+1. Hence we have a cover of Z* of size 
{t/2){£ + l)/2. Since t — [p— we have a cover with the properties we require 
and so the theorem follows in this case. 

Now suppose that £ is even. This implies that 2^/^ = —1 in Z*. Writing Z* as 
the union (4.1) of the cosets of (2), we see that 

C, = {ci2* I < i < £/2} U -{2*c, \ Q<t< £/2}. 

So, defining 

A = {22^c, |0< j< r^/41}, 
we see that Ct is covered by the [^/4] sets B2.p{d) with rf e Z?^. Defining 

£> = £>! U L>2 U • • • U A 
we have constructed a covering of Z* consisting of t\£/A\ sets B2^p{d). Since i 



(p — 1)/^, the theorem now follows by Theorem 3.1 □ 



Corollary 4.2. Let p be a prime such that p = 5 mod 8. Then 

f2{p)<l(^^+ip + 3)/'i. 

Proof. If p = 5 mod 8, then 2 is a quadratic non-residue modulo p. The quadratic 
residues are exactly those elements of Z* whose orders divide {p — l)/2, and hence 
4 must divide the order of 2. The corollary then follows by the final case of Theo- 
rem HITl □ 
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Corollary 4.3. The length f2{n) of the shortest n-ary 2-radius sequence satisfies 



2 V2, 

Proof. When n is sufSciently large, Theorem |2.3| imphes that there is a prime p 



2.1 



such that p = 5 mod 8 and n < p < n + rr-°^°. Since f2in) < f2{p) by Lemma 
the upper bound of the corollary follows by Corollary |4.2[ The lower bound follows 



by Lemma 2.2 □ 



5. Logarithmic functions, logarithms and tilings 

The results from Sections [s] and |4] indicate that good coverings (even better, 
tilings) of Z* by sets of the form B]. p{d) will give rise to short fc-radius sequences. 
This section defines the notion of a logarithm of length fc, and explains how to use 
these logarithms to construct certain high-dimcnsional tilings. In Section [6j we will 
use these tilings to construct good coverings of Z*, which in turn will produce short 
fc-radius sequences. 

5.1. Logarithmic functions and logarithms. Given A; > 1, a logarithmic func- 
tion (of length k) is a map / : {1, . . . , fc} such that f{ab) — f{a) + f(h) 
whenever a,b d {1, 2, . . . , fc} satisfy ab < k. A logarithm is a bijective logarithmic 
function. Logarithms are used to construct lattice tilings of n-dimensional space 
by semi-crosses; see Stein and Szabo [15^ Chapter 3]. They also arise in group 
theory, number theory and coding theory; see Galovich and Stein [S], Forcade and 
PoUington and Gordon |8|. 

Galovich and Stein [6] note that a logarithm of length fc exists whenever either 
fc~|- 1 is prime or 2fc-|- 1 is prime. To see this when fc-t- 1 is prime, let a be a primitive 
root modulo fc-|-l. Then the 'discrete logarithm' map, taking a E {1, 2, . . . , fc} to the 
element x G Zk such that a = mod fc -f 1, is a logarithm of length fc. When 2fc+l 
is prime, let a be a primitive root modulo 2fc + 1. Let g : {1, 2, . . . , fc} — > be the 
map taking a £ {1, 2, . . . , fc} to the element x G Z2k such that a = mod 2fc -I- 1. 
The composition of g with the natural surjection from Z2k to Z^, (sending x mod 2fe 
to X mod fc) is a logarithm of length fc. Computations by Galovich and Stein [B], 
and later by Forcade and Pollington [5^ , show that logarithms are common for small 
values of fc. In particular, their computations show that logarithms exist for all fc 
with fc < 204, except 195. We have verified and extended Forcade and Pollington's 
computations; we report our results in Section [9] 

Let p he a prime such that p = 1 mod fc, and let C be a primitive fcth root of 
unity in Z*. Let / be the map / ; {l,2,...,fc} — )• Zfc defined by the equation 
^{p-i)/k ^ (^fix) j^Qd p Then / is certainly a logarithmic function of length fc. If 
/ happens to be bijective, then / is a logarithm. We say (following Galovich and 
Stein [6]) that / is a Kummer-Mills- logarithm, or KM-logarithm, if it arises in this 
way for some prime p. (The Galovich-Stein definition for a KM-logarithm / is that 
/ extends to a fc-character Z* — >■ Zfc. The relevant fc-characters of Z*, restricted to 
{1, . . . , fc}, correspond precisely to the choices for our function f as ( varies over 
the primitive fcth roots of unity.) 

Theorem 5.1 (immediate from Theorem 3 of Mills [13]). Every logarithm of odd 
length k is a KM-logarithm. //fc = 2 mod 4 then a logarithm f of length k is a KM- 
logarithm if and only if f{m) is even whenever both m divides fc and m = 1 mod 4. 
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If k = mod A, then a logarithm f is a KM-logarithm if and only if f{m) is even 
whenever m divides k/A. If f is a KM-logarithm of length k, then there are infinitely 
many primes p that exhibit f as a KM-logarithm. 

A logarithm / of length k will be called a special KM-logarithm if either (i) k is 
odd, or (ii) k is even and /(m) is even whenever m divides fc/2. Any special KM- 
logarithm is a KM-logarithm, and indeed unless fc = 2, 4, or 6 mod 8, the concepts 
of a KM-logarithm and a special KM-logarithm coincide. (For k odd this is clear. 



If 8|fc and / is a KM-logarithm then from Theorem 5.1 we see that /(m) is even 
whenever m divides k/A. In particular this holds for prime values of to dividing k/A. 
Since 8|fc, the primes dividing fc/4 are the same as those dividing fc/2. Hence /(to) 
is even for each prime to dividing fc/2, and hence for all to dividing fc/2. Thus 
/ is a special KM-logarithm.) A special KM-logarithmic function will naturally 
be a logarithmic function that satisfies the parity conditions for being a special 
KM-logarithm, but need not be bijective. 

Lemma 5.2. There are infinitely many k for which there is a special KM-logarithm 
of length fc. Indeed we have: 

(i) if p = k-\-l is prime and 8|fc, then there is a special KM-logarithm of length fc; 

(ii) if p = 2fc -|- 1 is prime, then there is a special KM-logarithm of length fc. 

Proof, (i) For such p any KM-logarithm is a special KM-logarithm, and a discrete 
logarithm in Z* is a KM-logarithm. 

(ii) Let / be a logarithm of length fc induced by a discrete logarithm in Z*. 
If fc is odd, then we are done: / is a special KM-logarithm. If fc is even, then 
p = I mod 4, and hence for any odd prime q dividing fc we have = = 1, 
and since fc is even we conclude that f{q) must be even. If 4 | fc then p = 1 mod 8, 
so — 1, and hence /(2) must be even. We see that in all cases / is a special 
KM-logarithm. □ 

5.2. Logarithms and tilings. Let r be a positive integer, and let C C Z*" be a 
finite set; we refer to C as a cluster. A Z-tiling of 7/" by the cluster C is a set L <ZJ7 
of points such that the union UxslI^ + ^) ^ disjoint union and is equal to 17 . If 
L is an integer lattice of determinant fc (i.e., L is a subgroup of U of index fc), then 
L is a Z-tiling of T7 by C if and only if C is a set of coset representatives for L\ in 
other words if and only if the natural map from U to the quotient W jL is bijective 
when restricted to C. The tiling terminology comes from imagining C as a set of 
r-dimensional unit hypercubes, indexed by the positions in TH" of their centres. A 
Z-tiling is then a tiling of W!" by translates of C in the usual sense. 

We now show how to use logarithms to produce certain tilings. Let qi,q2, . ■ . 
be the primes in the set {1,2,..., fc}; so r = 7r(fc). Let No be the set of non-negative 
integers, and let Ck be the cluster of points in Z*" defined by 

(5.1) = |(zi,*2,...,z.)eN5 I H'/j' <fc 

We identify the points in Ck with the integers 1, . . . , fc. 

Theorem 5.3. Let k be a fixed positive integer, let r = 7r(fc) and define Ck as 
above. If there exists a logarithm of length k, then there exists a 'Z-tiling of IT by 
Ck- 



8 



SIMON R. BLACKBURN AND JAMES F. MCKEE 



Proof. Let / be a logarithm of length k. Define a group homomorphism : Z*" — > 

by 

(^((Zl, ^2, . . . , ir)) = ilfiqi) + i2.f{q2) + h irfiqr), 

so (f> sends the jth standard basis vector to f{qj)- 

Define the lattice L C to be the kernel of 4>. By our remarks above, the union 
Uajsil^ + Cfc) is a tiling if and only if is a bijection when restricted to Ck- Now, 

^((n, 12, . . . , ir)) = fiq^ql" ■ ■ ■ q'/) 

for all (ii, ?2, ■ ■ ■ ,ir) G Ck, since f{ab) = f{a) + f{b) whenever a, b, ab £ {1,2,..., fc}. 
So the restriction of cj) to Ck is equal to /. Since / is a bijection, we have the tiling 
we require. □ 



6. /c-RADIUS SEQUENCES 

This section contains the main result of our paper: the leading term of fk{n) for 
infinitely many values of k (including all k < 194). We need to establish a simple 
geometric lemma before we begin the proof. 

For linearly independent vectors xi,X2, ■ ■ ■ ,Xr G M'', we define the r-dimensional 
parallelotope P{xi, X2, ■ ■ ■ , x^) to be the region 

P{xi,X2, ... ,Xr) = l^ttiXi I < flj < l| . 

Lemma 6.1. Let r be a fixed positive integer, and let S be a positive real number. 
For any linearly independent vectors Xi,X2, . . . ,Xr € K'' with \xi\ < \x2\ < ■ • • < 
\xr\, let P be the parallelotope P ~ P(xi,X2, ... ,Xr), and let Q be the set of points 
in W within a distance 6 of P. Then 

Vol,(Q) < Vol,(P) + l^fjl^^l) ' 

where the implicit constant depends only on r and S. 

Proof. The parallelotope has a bounded number of proper subfaces F (it has (^)2''~' 
subfaces of dimension i). For a proper subface F, define Qp to be the set of points 
w e M'' such that (i) v is within a distance S of F, and (ii) the nearest point in F 
to V does not lie in a proper subface of F. Then 

Q = Pu[jQF, 

F 

where the union runs over all proper subfaces F of P. This decomposition in the 
case r = 2 is illustrated in Figure [l] Here, the parallelotope P (a parallelogram as 
we are in dimension 2) is in bold; the circles are the sets Qp when F has dimension 
and the rectangles are Qp when F has dimension 1. This decomposition shows 
that 

(6.1) Vol^(Q) < Vol^(P) +^Vol^(gj.). 

F 

It is not difficult to see that when F is i-dimensional, Qp is congruent to FxBj.-i, 
where B,.-i is the (r — «)-dimensional ball of radius S. Now F is congruent to an 
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Figure 1. Decomposing Q 



i-dimensional parallelotope, spanned by a proper subset S C {1,2, ... ,r} of the 
vectors x,-. So 



voi,(F) < n < n i^.i<ni^^- 

Since the baU Br-i is contained in a hypercube with sides of length 2(5, we find that 
Yolr-iiBr-i) < i26y. Hence 

r I 

Vo\r{QF) = Vol,(F) X Vol,_,(B,_0 < {25y [] \x,\ = O [] |.t 

J=2 \j=2 



The lemma now follows by (6.1), since the number of terms in the sum in (6.11 is 



bounded. □ 
Theorem 6.2. Let k be a fixed positive integer, and let r = 7r(fc). Suppose there 



exists a 'L-tiling L of U by Ck (as defined by (5.1 1). Then 



I (n\ „, 2 



fk{n) = -^^j+0{n'/\ogn). 

{Recall that such a tiling exists whenever there is a logarithm of length k, so in 
particular when k < 194 or when either k + 1 or 2k + 1 is prime.) 



Proof. By (1.1) we may assume that k > 1. Let qi,q2,...,qr be the primes in 
the set {1, 2, . . . ; k}, where qi — 2. Let p be any prime such that (^^^ — ^1 

( i^) ~ ^ ^'^^ alH e {1, 2, . . . , r}. A sufhcient condition for a prime p to satisfy these 
character constraints is that p = —1 mod 85293 ■ ■ ■ q^ so by Theorem |2.3| we find 
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that for all sufficiently large n there exists such a prime p with n < p < n+ [12 \ . 
We aim to show that 

Mp)<l(^^+o{py\ogp). 

This suffices to prove the theorem, since we may then deduce that 

Mn) < Mp) < ^("+ +0((n + nO-^2^)Vlogn) 

404-O(nVlogn), 

as required. 

Recall the definition of the sets Bk,p{d) from Section [sj We will establish the 
bound we require by constructing a good covering of Z* by sets of the form Bk^p{d) 
and then applying Theorem |3.1[ 

Define the subgroup H of the multiplicative group Z* by H = {qi, q2, . . . , qr)- 
Define £ — \H\. We may write Z* as the disjoint union of the cosets of H: 

z; = CiuC2u---uCt, 



4.1 



where Ci — CiH for some elements Ci € Z*. Just as in the proof of Theorem 
we note that | = ^ for alH G {1, 2, . . . , and f = (p - 1) /i. 

Since ^^i^ = 1 for i e {1, 2, . . . , r}, we see that H consists of quadratic residues 

modulo p and so all the elements in a coset share the same Legendre symbol. For 

any i G {1,2,..., t}, —d — Cj for some j G {1,2,..., t}. Since j = — Ij the 

Legendre symbols of the elements of Ci have opposite sign to the Legendre symbols 
in Cj, and so in particular i ^ j. Hence t is even and, by relabelling cosets if 
necessary, we may assume that 6*2^+1 = —C2S for all s G {1,2,..., t/2}. 
For d G Z;, define 

Ak.,p{d) - d{l, 2, . . . , fc} = {d, 2d, . . . , kd}. 

Note that Bk.p{d) — Ak.p{d) U —Ak,p{d). We aim to show that there exists an 
integer w such that w < {£/k)(l + 0(l/logp)) with the following property: there 
exist elements di, o?2, ■ • • , d™ € H such that the sets Ak.p{di) cover H. This is 
sufficient to prove the theorem; we can see this as follows. Since Bk,p{d) — Ak,p{d)U 
—Ak,p{d), the w sets Bk,p{di) cover H U —H. So the wt/2 sets Bk.p{c2sdi) where 
i G {1, 2, . . . , fc} and s G {1, 2, . . . , t/2} cover Z*. Theorem 3.1 now implies that 



fin 

(p — l)w , , , 

_(p-i)(i + o(Vi"gp))(pifc) 



k\2 



2k 

O(pVlogp), 



as required. 

It remains to show the existence of the elements di , d2 , . . . , d^, G H with the 
properties we need. Let L C Z'" be such that {z + : z G L} is a Z-tiling of 1/ . 
Let (/) : Z'' — > be the group homomorphism sending the ith standard basis vector 
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of Z'' to Qi for «e{l,2,...,r}. Note that is written additively but H is written 
multiplicatively, so (j){x + y) — 4){x)4){y) for all x,y G . Also note that, using the 
definition of Cfc, 

(j>{x + Ck) = Ak,p{(l){x)). 
Let K = kcrt/), so i^T is a lattice. Since H is generated by qi,q2,. . . ,qr, we see 
that (p is onto and so K has determinant £. Let Xi, X2, ■ ■ ■ , Xj. be a basis for this 
lattice. We may choose our basis such that 

r 

n l^^'l ^ 2'-('-i)/4det(i^) = 2''('-i)/'*^, 
1=1 

by LLL- reducing this basis if necessary; see Cohen |31 Section 2.6], for example. 
Moreover, by reordering the basis if needed, we may assume that \xi\ < \x2\ < 

• • • < \Xr\- 

We claim that the shortest non-zero vector in K has length at least (logj,p)/r, 
and so in particular |xi| > {\ogf.p)/r. To establish our claim, suppose that the 
vector (zi, i2, . . . , v) G K has length less than (logj,p)/r. We need to show that this 
vector is zero. Our bound on the length of the vector implies that < (logj.p)/r 
for all j, and so J2]=i Njl < logfcP- Define non- negative integers i'j and i" by 
ij — max(ij, 0), i" = — min(ij, 0). Then 

(ii,i2, • ■ • ,«r) - («i,«2, ■•■,«") = (*i,«2, . . . , v) e ^ 

and so 

(6.2) qfqf ■■■ir^ i^, • • ■ , z^)) = z^', ■ • ■ , i';)) = qUi ' ' ' <h^ mod p. 

But both the left hand side and right hand side of the equation above are integers 
bounded above by p. For example. 



So no modular reduction can have taken place in (6.2 1 and therefore i'^ = i" for 



j € {1, 2, . . . , r}. This implies that ii = 12 = ■ • • = v 0, and our claim follows. 
Let it! C Z'' be defined by 



{y^T- \ y = Y a,x, : < a, < 1 



Since R is the intersection of 1/" with a fundamental region of the lattice we 
see that = H (indeed, the restriction of (/) to i? is bijective). Define W to be 

the set of translations z G L in our tiling such that {z + Ck) r\ R ^ %, and define 
w = liyl. We have that 



U ((/.(z + O) = 0( U {z + Cu)] 



D 0(i?) = H, 



and since (j>{z + Ck) = Ak^p{(l>{z)) we have a covering of H consisting of w sets of 
the form Ak.p{d). To prove the theorem, it suffices to show that w < {£/k)(l + 

o{i/iogp)y 

Define Ck to be the union of those unit hypercubes in Z*" whose centres lie in Ck ■ 
Let 6 be the diameter of Ck', note that 5 depends only on k. Now, R is contained in 
the parallelotope P — P{xi,X2, ■ ■ ■ , x,.) of volume £. Define Q C M'' to be the set 
of points within a distance 6 of P. When z e W, the set z + Ck contains a point of 
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R and so 2 + Cfc C Q. Since W C L and L is a Z-tiling of Ck, the interiors of the 
sets z +Ck with z £ W are disjoint, and so Yo\r{Q) > J2zew^'^^r{x + Ck) — wk. 
However, Lemma |6. 1 1 imphes that 



Vol^(g) < Vol^(P) + [W 



Xj 



\i=2 




the last line following since nl=i 1^*1 = 0{t) and > (logj,p)/r. Hence 

^voi,(g)< (^/fc) (1 + 0(1/ logp)), 

as required. □ 

7. fc-RADIUS PRIMES 

We now introduce the notion of a A:-radius prime, and show that fc-radius primes 
exist if and only if there exists a special KM-logarithm of length k. This allows us 
to prove (Corollary 
that of Theorem 16.2 



7.3 1 an upper bound on fk{n) with a better error term than 
provided we are willing to accept a bound for infinitely many 
values of n rather than for all values. In ^ we will determine the density of fc- 
radius primes, and this will allow us, for values of k where fc-radius primes exist. 



to improve the error term in Theorem 6.2 without restricting n; see Theorem 



We say that a prime p is a k-radius prime if the following two conditions both 
hold: 

(7.1a) • p=\ mod 2fc ; 

(7.1b) • the elements l^P-i)/'^', 2(^-1)/^ . . . , fcC^-i)/*^ in Z; are pairwise distinct. 
Note that is a fc-th root of unity in Z* for any 2; e Z*. So the condition 



(7.1b I is quite a strong one, as there are only fc such roots of unity. We remark, for 



later use, that the congruence condition (7.1a) implies that {p — l)/fc is even. 



Proposition 7.1. Let p be a k-radius prime. Then there exists a p-ary k-radius 
sequence of length ((p — l)/2fc)(p + fc — 1) + 1. In particular, 

fkiP)^u^^+oip) 



k V2y 

when p is a k-radius prime. 

Proof. Let a be a primitive root modulo p. Define D C Z* by 

D = {a''^ :0<i< (p- l)/2fc}. 
We claim that the sets Bk,p{d) where d G D cover Z* Since D C Z* has size 



(p — l)/2fc the proposition will then follow by Theorem 3.1 It is sufficient to show 
that our sets are disjoint, for then their union has size 2fc|Z)| = p — 1. Suppose, for 
a contradiction, that there exist < ii < 12 < (p — l)/2fc such that Bk,p{a'^^^) n 
Bk,p{ct-^^^) 7^ 0. By the definition of Bk,p{d), there exist ri, r2 € {1,2,..., fc} such 
that ria''*! = ±r2a'^*2 mod p. Raising both sides of this equation to the power 
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{p ~ gives the equation ^^^^ = (il)'?" '^)liij.'yV 1)/'= jjiodp. Since p is a fc- 
radius prime, the congruence condition ( 7.1a[ ) imph es tha t (p — l)//c is even. Hence 



rf ^^^^ = ^'^'^ mod p. But now the condition (7. lb I imphes that ri = r2, and 



hence that a'^'^ = ±a'^'^ mod p. If a''^^ = a'^'^ modp then clearly ii = 12 and we 
have our contradiction. If a^^^^ = — a*^*^ mod p then 

kii — A:i2 ^ {p ^ mod (p — 1). 

Again we have our desired contradiction, since our bounds on ii and 12 imply that 
k{ii — 12) has absolute value less than (j> — l)/2. □ 

The next theorem provides a connection between fc-radius primes and special 
KM-logarithms. It shows us that for infinitely many values of fc there exist fc- 
radius primes. In particular, if 2fc + 1 is prime, or if fc is divisible by 8 and fc + 1 is 
prime, then there exist infinitely many /c-radius primes. 

Theorem 7.2. I] there is a special KM-logarithm of length k, then there are infin- 
itely many k-radius primes. Conversely, if there is a k-radius prime then there is 
a special KM-logarithm of length k. 

Proof. Let / : {1, . . . , fc} ^ Zfe be a special KM-logarithm. 



If k is odd, then Theorem 5.1 gives infinitely many primes p = 1 mod k satisfying 



the condition ( 7.1b[ ), and with k odd these automatically satisfy the congruence 



condition (7.1a) too. 

If k is even, then define / : {1, . . . , 2k} — >■ Z2fc to be any logarithmic function that 
agrees with / modulo k on the primes between 2 and k. In particular, we then have 
that /(a) = /(a) mod fc for 1 < a < fc. Since / is a special KM-l ogar ithm, the values 



of / provide legitimate character values mod 2k for Theorem 5.1 to apply: there 



exist infinitely many primes p = 1 (mod 2k) such that a^^ i)/2fc ^ ^fi<^) u^odp for 
1 < a < k, where Z2fc is some primitive 2fc-th root of unity in Z* We certainly 



have the congruence condition (7. la I, but now we must check the condition (7.1b). 
Suppose that o^p-i)/*-' = b'-P-^^^^od p. Then a'-P'^^^'' = ±b^P-^^/'^k j^odp, so 

z^Jj^^ = 'i^z^J^^ mod p. Since —1 = z^p. modp, we deduce that f(a) = f{b) mod k. 
But / agrees with / mod fc, and / is a logarithm, so we must have a ^ b. 

Conversely, suppose that p is a fc-radius prime, and let Zk be a primitive k-th 
root of unity in Z*p. Define / : {1, . . . ,A:}_> by a^P-^)/*^ = z^^"^ mod p. Then 



/ is certainly a logarithm. We now use (7.1a) to proceed via quadratic reciprocity 



(as in the proof of Lemma 5.2) to verify that / is a special KM-logarithm. 

□ 

Corollary 7.3. Let k be a positive integer. Suppose there exists a special KM- 
logarithm of length k. (Recall that this assumption holds when 2fc -t- 1 is prime, or 
when k + 1 is prime and 8\k.) Then fk{n) — \ +0{n) for infinitely many values 
ofn. 

8. The density of fc-RADius primes 

We believe the notion of a fc-radius prime is a very natural one. So it is of interest 
to determine the density of such primes; and this can be done, using some of the 
theory of cyclotomic field extensions. For technical background material, see |10j . 
As an application of our result on the density of fc-radius primes, we improve the 
error term in Theorem |6. 2 1 for values of fc such that special KM-logarithms exist. 
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Theorem 8.1. Let k be a fixed positive integer, and let fspec{k) be the number of 
special KM-logarithms of length k. Let 



(8.1) 



1 



/spec (^) 

Lp{2k) ' fc^(fc) 



Cfe = < 



1 



/,pec(fc)2"('=/2) 



if k is odd. 



if k is even. 



[.Lp{2k) k^f") 

There exists a positive constant Ak such that the number of k -radius primes less 
than or equal to x is 

X 



Ck 



logo; 



O (xcxp{-Ak\/\ogx) 



as X ^ oo, where the implied constant depends on k. 

Heuristically the main term is clear: the fraction l/(p(2k) gives the proportion 



of primes satisfying (7.1a); the remaining fraction gives the proportion of special 
KM-logarithmic functions that are special KM-logarithms (when k is even, the 
parity constraint on each prime dividing k/2 means that the number of special 
KM-logarithmic functions is only k^i^) /2'^'^^/'^^ rather than k'^^^^). We shall make 
use of a theorem of Elliott (H Theorem 1] . See also [H] for a discussion of similar 
density results, some of these conditional on the generalized Riemann Hypothesis. 

Proof. Suppose first that k is odd. Let C,k — exp(27ri/fc). We work in the cyclotomic 
field Q(Cfe), and associate any logarithmic function / : {1, . . . , fc} — with the 

vector e, = (Cf\Cfe^(^',...,Cr^ 

If p = 1 mod fc, then p splits completely in Q(Cfe) as a product of ip{k) degree-1 
prime ideals. Apart from a finite number of ramified primes, all degree-1 prime 
ideals in Q(Cfe) arise in this way. For any degree-d prime ideal p, the fc-th power 
residue symbol ^ is defined for j3 e 1i[Ck] to be zero if p | /3, and otherwise 

is the unique power of C,k that is congruent to /J^p''^^)/'^ mod p. We can associate 
to any prime ideal p the vector £p = ^ • Let S be the finite set 

of all vectors that arise as Ep for some degree-1 prime ideal p dividing a rational 
unramified prime p. For a degree-1 prime ideal p dividing p = 1 mod fc, we observe 
that the powers of Cfe appearing in Cp correspond to the powers of some primitive 



fcth root of unity modulo p that appear in the sequence in (7.1b I. The different 



prime ideals dividing p correspond to different choices of the fcth root of unity. For 
condition (7.1b) to hold we require the entries of £p to be distinct. 

Given any vector e = (ei, . . . , e^) whose entries are powers of Cfe, Elliott's theorem 
[31 Theorem 1] gives an asymptotic estimate for the number of prime ideals p with 
norm less than or equal to x such that £p = e, namely 

N{k,k) 



k^ 



-Tr{x) + O ^xexp(— yli/loga;)^ , 



for some positive constant A, where 
(8.2) 7V(A:,fc)=^(£r 
and where the sum in 



r 



2 ) is over all 1 < i^i, . . . , i/fe < fc such that ■ ■ ■ k"" = /?" 



for some /S e I'lCk]- Since the number of prime ideals of degree greater than 1 and 



CONSTRUCTING fc-RADIUS SEQUENCES 



15 



k 


1 


2 


3 


4 


5 


predicted 
observed 


1.00 
1.00 


0.250 
0.250 


0.111 
0.111 


0.00 
0.00 


0.00160 
0.00161 


k 


6 


7 


8 


9 


10 


predicted 
observed 


0.00463 
0.00464 


0.00250 
0.00250 


0.000977 
0.000974 


0.000610 
0.000600 


0.000200 
0.000202 



Table 1. Predicted and observed (primes up to 10^) densities of 
/c-radius primes 



norm less than or equal to x is dominated by the error term in this estimate, we 
have the same asymptotic formula if we restrict to degree-1 prime ideals. Now we 
observe that ii e ~ Sp is in S, then for each term in (8.2) we have 



Hence each e G S gives us the same main term in the asymptotic estimate for the 
number of degree-1 prime ideals of norm less than or equal to x such that £p = e, 
and with the same form for the error term. Now Mills' thoerem [131 Theorem 3] tells 
us that the e £ S are precisely those that correspond to (special KM-) logarithmic 
functions (as we are assuming k is odd there are no parity constraints), hence 
amongst all degree-1 prime ideals the proportion that yield special KM-logarithms 
(corresponding to fc-radius primes) is precisely fspec{k)/k^^''\ 

For fc even there are a few twists, but essentially the same argument works. We 
now work in Q(C2fc)- For each special KM-logarithmic function / : {1, . . . , fc} — )■ 
we associate a set of 2^^*^^ vectors (ei, . . . , e^), where each Si is a power of C2fc such 
that for primes q between 1 and fc we have either Sq = c/i'^'' — Cii^''^^ ~ 

—C2k'^\ extended logarithmically to all other £j. The parity conditions for / being 
a special KM-logarithmic function are precisely those needed for Mills' theorem 
to tell us that for any of these 2'^^*''^ vectors e there exist infinitely many degree-1 
prime ideals p € Q(C2fc) such that ^ = ^p = ((p) '■■■'(?) )■ Bearing in mind 

these parity conditions, we have 2'^*^'^-'fc'^'-'^''/2"'-*''/^-' vectors associated to special 
KM-logarithmic functions, and the proportion of degree-1 prime ideals that corre- 
spond to special KM-logarithms is /spoc(fc)2"'-'^^^V^'^'"'^''- The rest of the argument 
proceeds as before, using N{2k, fc) rather than (fc, fc) (the vj run between 1 and 
2fc). □ 

Corollary 8.2. The proportion of all primes less than or equal to x that are k- 
radius primes is c^, defined by (18. 1|). 



To illustrate this corollary. Table [T] compares the predicted densities with the 
proportion of fc-radius primes observed for primes up to 10^ and fc < 10. 



Corollary 8.3. If fspedk) > 0, then with as in Theorem 8.1 there is a k-radius 
prime between n and n + O (n log n exp(— yl^ \/\og n)) . 
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Proof. From Theorem |8.1| we have that the number of /c-radius primes between n 
and n + g (with 1 < g < n) is 

\T~r~^~^ ^ r~] + 0{nexp{-Aky/\ogn)) . 
\\og{n + g) \ognJ 

This is strictly positive provided that the main term Ck iog"^+g) — Togn) strictly 
larger than the absolute value of the error term. This main term is asymp- 
totic to Cfeg/logn as n — oo, so with g larger than some constant multiple of 
n log n exp(— Afc-^/log n) we are assured of a fc-radius prime between n and n+g. □ 

Theorem 8.4. Let k be a fixed positive integer such that a special KM-logarithm 
exists. {Recall that a special KM-logarithm exists, in particular, when 2fc + 1 is 



prime.) Define the positive constant as in Theorem 8.1 Th- 



fk{n) = +0{n^ log n exp (- Afe ^^log n)). 



We remark that this improves the bound of Theorem 6.2 for those values of k 



where a special KM-logarithm exists, since exp(Afc-v^log n) grows faster than (log 
Proof. Let p be the smallest fc-radius prime such that n < p. Then p = n 



where g = O (n log n exp{~Aky/\og n)) , by Corollary 8.3 Proposition 7.1 implies 
that /fc(p) = i(P)+0(p) and so 

fkin) < fu{p) = ^ (2) + 0{P) - ^ (2) + 0{ng), 
as required. □ 

9. Computational results 

In this section, we present our computational results. We begin by describing the 
techniques we used to design an efficient search for logarithms, and to enumerate 
them. We then describe the results themselves. 

9.1. Finding logarithms. Let gi < • • • < g,. be the primes in {l,...,fc}. A 

logarithmic function / of length k is uniquely determined by the values of /(?i), 
. . . , f{qr), which can be chosen independently, so there are fc'^^*^^ — k^ of them. 

Let V = 'L}^ be the set of all vectors of length fc, with components in Zfc. We 
can conveniently represent any logarithmic function / of length k as the element 
V/ = (/(I), /(2), . . . , /(fc)) e V . For each prime < k, define e^ € F by letting 
the jth component of e, be the exponent of the largest power of g^ dividing j. 
Then / is a logarithmic function with /(gi) — ai for i G {1, 2, . . . , r} if and only if 
V/ = X]i=i o.i^i- So the set of logarithmic functions is a Z^-submodule of V , with 
generating set ei, e2, . . . , e^. In what follows, we identify a logarithmic function / 
with its representation vj. 

Logarithms are bijective logarithmic functions, so a naive strategy would be to 
search over all A;'" logarithmic functions. Forcade and PoUington [S] observed that 
it is sufficient to search over all possibilities for /(gi) with qi < k/2. Once these 
values are fixed, / is determined on all (fc/2)-smooth numbers in {l,2,...,fc}. 
If the resulting partial function is injective with image X, then / may always 
be completed to a logarithm by assigning the remaining values of /(g^) for g^ > 
k/2 to be an arbitrary permutation of Z/j \ X. This trick reduces the size of our 
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search from k'"^''^ to k^'^'^^^h Another way of thinking of this is as follows. We 
partition the set of logarithmic functions into parts; two functions are in the same 
part if and only if they agree on qi for Qi < k/2. We search over parts rather 
than logarithmic functions. Each part contains (7r(fc) — 7r(fc/2))! logarithms if the 
appropriate partial function is injective, otherwise the part contains no logarithms. 
In our computations, we made use of a generalisation of this, which we now describe. 

Partition the primes gi , . . . , q,. into Mocks Bi, . . . , Bf in following way. Any 
prime < y/k appears on its own in a block of size 1. Primes q in the interval {Vk, k] 
are grouped together in blocks according to the value of i such that fc/(^+l)<g< 
k/£. Observe that if tt is a permutation of {1, ... , r} that induces a permutation of 
each block, and / — X]i=i '^i^i is a logarithm, then so is X]i=i Logarithms 
that are related in this way we call combinatorially equivalent. 

If / is a logarithm of length k, then so is af for any a G . Two logarithms / and 
g will be called linearly equivalent il f = ag for some a G Z^. Two logarithms / and 
g will be deemed equivalent if there is a logarithm h such that / is combinatorially 
equivalent to h and h is linearly equivalent to g. For finding and counting logarithms 
it is extremely useful to work up to equivalence. 

We remark that a small modification of this technique works when counting 
KM- logarithms or special KM-logarithms. If / is a KM-logarithm (or a special 
KM-logaritlini) , then so is af for any a G Z^, for if k is odd then there arc no 
parity conditions to check, whilst if k is even then a must be odd, and the parity 
of a/(j) (reduced modulo even k) is the same as the parity of /(j). When dealing 
with special KM-logarithms or KM-logarithms, any primes dividing k must be put 
in singleton blocks. Since for simple logarithms the only primes q in blocks of size 
greater than 1 satisfy q > \fk, there is at most one block that is affected by this 
change, having one of its primes removed to a singleton block. 

Given fc, our search for a logarithm / proceeds by attempting to assign values 
for f{qi), . . . , f{qr) in that order. Having found compatible values for f{qi), ■ ■ ■ , 
f{qj-i) and having computed Vj_i = f{qi)ei, we consider each possibility 

for f{qj) in turn: the vector Vj = Vj-i + f{qj)ej is computed; the gj-smooth 
components of Vj are extracted (using a precomputed list of gj-smooth numbers), 
sorted, and checked for repeats; if a repeat is found then the next value of f{qj) is 
tried, backtracking once all possibilities are exhausted; if no repeat is found then j 
is increased. 

The backtracking makes use of equivalence. Within each block {qi^ ,qi^} of 
primes, we constrain the values of the logarithms to satisfy f{qii) < ■ ■ ■ < /(ftj. 
Linear equivalence is exploited by insisting that /(2) | k, and restricting values of 
/(3) so that there does not exist a with gcd(a, k) = 1 and af{2) = /(2) mod k such 
that o/(3) mod k is smaller than /(3). 

For computing KM-logarithms and special KM-logarithms, the parity conditions 
are checked on g^-smooth numbers before accepting the legitimacy of f{qj). 

9.2. Counting logarithms. To count the total number of logarithms (or KM- 
logarithms, or special KM-logarithms) of length fc, we need to know the size of each 
equivalence class, and how many representatives of each class appear in the output 
of our algorithm (now not terminating when a logarithm is found, but continuing 
to search the entire space of logarithmic functions). If / is a logarithm, scaled so 
that /(2) I k, then the number of a S Z|. such that a/(2) = /(2) mod k equals 
the number of a G with a = 1 mod {k/f{2)), which is the size of the kernel of 
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n 


1 


2 


3 


4 


5 


6 


7 


/log(fc) 


1 


1 


2 


2 


8 


10 


36 


,/spcc (^) 


1 


1 


2 





8 


2 


36 


n 


8 


9 


10 


11 


12 


13 


14 


fiog{k) 


40 


24 


20 


140 


136 


936 


624 


/spec (^) 


16 


24 


8 


140 





936 


312 


n 


15 


16 


17 


18 


19 


20 


21 


/log(fc) 


416 


256 


3 648 


2 088 


30 240 


16 704 


9 792 


/spec (^) 


416 


96 


3 648 


576 


30 240 


4 608 


9 792 


n 


22 


23 


24 


25 


26 


27 


28 


/log(fc) 


9 000 


103 488 


86 400 


72 960 


36 576 


22 896 


12 096 


/spec(fc) 


3 360 


103 488 


10 368 


72 960 


13 752 


22 896 


1,296 


n 


29 


30 


31 


32 


33 


34 


35 


/log(fc) 


134 400 


105 216 


2 671200 


1708 800 


794 400 


396 288 


145 152 


/spec(fc) 


134 400 


23 424 


2 671200 


556 800 


794 400 


202 752 


145 152 


n 


36 


37 


38 


39 


40 


41 


42 


/log(fe) 


109 440 


3 594 240 


2 244 672 


1202 688 


660 480 


17606 400 


16104 960 


/spec (^) 


7488 


3 594 240 


1013 472 


1202 688 


102 912 


17606 400 


2 021760 



Table 2. The number of logarithms and special KM- logarithms 
of length k 



the natural homomorphism Z*^ ^fe//(2)' namely f{k)/tf[k/f{2)'^. To ease the 
counting, we remove the constraint on /(3), and check for each logarithm that 
the vector of logarithm values for the primes in singleton blocks scales to precisely 
ip{k)/ip{k/ f{2)') distinct vectors on multiplication by elements of Z^. If the primes 
below k fall into blocks of lengths £i , . . . , , then each logarithm / produced by 
our algorithm contributes ip[k/f{2)) JlJ^i ^i' to the total number of logarithms. 

9.3. Results and comments. We write fiog{k) for the number of logarithms of 
length fc, and /spcc(fc) for the number of special KM-logarithms. Table [2] shows 
/log(fc) and /spcc(fc) for small values of k. 

Table [2] shows that logarithms of length k are common when k is small. This 
entire table took just a few minutes to compute. However, logarithms of length k do 
not always exist: Forcade and PoUington [5] compute that there are no logarithms 
of length 195, but that logarithms of length k exist for all other values of k with 
k < 204. Their computations also show that there exists no logarithm of length 
255. 

We pushed the computations for finding logarithms, KM-logarithms, and special 
KM-logarithms up to A; = 300. There is no logarithm of length k ioi k = 195, 205, 
208, 211, 212, 214, 217, 218, 220, 227, 229, 235, 242, 244, 246, 247, 248, 252, 253, 
255, 257, 258, 259, 263, 264, 265, 266, 267, 269, 271, 274, 275, 279, 283, 286, 287, 
289, 290, 291, 294, 295, 297, 298. For all other k < 300 there exists a logarithm of 
length k. In addition to the above list, there is no KM-logarithm of length k for 
k = 184, 234, 236. In addition to all of these, there is no special KM-logarithm for 
fc = 4, 12, 60, 180, 182, 190, 196, 222, 238, 268, 276, 282, 292. When no logarithm 
existed, the computations for larger values of k typically took about 50 hours each. 
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p 


2 


3 


5 


7 


11 


13 


17 


19 


23 


29 


31 


37 




1 


122 


29 


100 


48 


108 


113 


239 


190 


166 


176 


64 




p 


41 


43 


47 


53 


59 


61 


67 


71 


73 


79 


83 


89 


f{p) 


18 


134 


145 


226 


38 


72 


81 


43 


182 


196 


249 


253 



P 


97 


101 


103 


107 


109 


113 


127 


131 


137 


f{p) 


23 


54 


120 


184 


217 


233 


260 


262 


271 



Table 3. A logarithm of length 277 



using PARI/GP [13] on a laptop, with the worst case by far being k = 294 (which 
took about 150 hours). 

Evidence from the enumeration in Table [2j together with heuristic combinatorial 
considerations, seem to point towards logarithms of length k being common when 
k is prime. However, this is not the case in general: the smallest prime k for which 
there is no logarithm is k = 211. 

The largest known k for which there exists a logarithm without either k + I 
or 2/c + 1 being prime is fc 277. For k above 200, these 'sporadic' logarithms 
appear to be increasingly rare: k = 201, 202, 203, 206, 207, 213, 223, 225, 234, 
236, 237, 241, 272, 277. (Could it possibly be the case that the only large values 
of k such that logarithms exist are when fc + 1 or 2fc + 1 is prime?) To illustrate 
our output, we present the values of a logarithm of length 277: it is enough to give 
the values at primes below k/2 (which determine the values at all (fc/2)-smooth 
numbers up to k), as the values at remaining primes below k may be assigned the 
missing logarithm values in any order. The bold numbers in Table |3] indicate the 
start of each of the blocks of primes; within each block, the logarithm values are 
seen to increase. 

Forcade and PoUington remarked in [5 that whenever they had found a logarithm 
/ of length A:, they had managed to find one with /(2) = 1. In fact there are many 
exceptions to this. For k < 288, the exceptional values are k = 176, 198, 200, 215, 
216, 219, 224, 225, 232, 234, 240, 249, 250, 256, 260, 262, 276, 280, 282: for these 
values of k there exist logarithms of length fc, but none with /(2) — 1. 

10. Open problems 

The big open questior{^for the asymptotics of fc-radius sequences is the following: 

Open Problem 1. For any fixed k, determine lim„_i.oo /fe ('t-) / (2) ("if exists). 

The results in this paper show that 1/fc < fkH/i'^) < l/(fc + o(fc)) for suffi- 
ciently large n, but is it the case that lim„_j.oo /fe («) / (2) — In particular, is 
it the case that lim„_j.oo /195 («)/ ( 2) = 1/195? Our results say that this limit (if 
it exists) lies in the interval [1/195, 1/194]. It is not difficult to see that whenever 
there exists a covering of Z'^'^^^^ by translated copies of Cigs with density 1 + e. 



Recent work by the first author has used probabiUstic methods to prove that this limit exists 
and equals However the problem of explicitly constructing asymptotically good fc-radius 

sequences for all fc remains open. 
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then the hmit is at most (1 + e)/195. So one approach to determinmg this hmit 
might be to find good coverings of Z'^^^^^) by copies of C195. 

Open Problem 2. Determine the set of integers k such that there exists a loga- 
rithm of length k. 

Is it the case that for all sufficiently large k, logarithms exist if and only if 
either A; + 1 or 2fc + 1 is prime? If this is not true, is it the case that the set of 
integers k such that there exists a logarithm of length k has density zero? Our 
computer searches for logarithms suggest this is very possible. In particular, is 
there an infinite sequence of integers k such that no logarithms of length k exist? 
We believe this is very likely, and give a heuristic argument below. 

Let Mfc be the largest size of the image of a logarithmic function from {1, . . . , fc} 
to Zfc, and let Rk < k he maximal such that there exists a logarithmic function 
from {1, . . . , fc} to Zfc whose restriction to i?j.-smooth elements is injective. Thus a 
logarithm exists if and and only if Mk — Rk = k, and when a logarithm does not 
exist the numbers Mk and Rk give natural measures of failure. 

Open Problem 3. Investigate Mk and Rk- 

We have that lim supj^^^ Mk/k — lim supj^_^^ Rk/k — 1. Is it the case that 
liminffc_j.oo -^fe/fc = liminffe_j.oo ^fe/fc = 0? We give a heuristic argument for Rk 
(with thanks to Carl Pomerance for suggesting this line of reasoning). Fix small 
e > 0. Let y) denote the number of y-smooth numbers below x. With u = 
log a;/ log 2/ we shall use repeatedly the estimate 'ip{x,y) — a;u^"+°^"\ uniformly as 
It — > 00 with u < (1 — e) log a:/ log log X [T]. Let the rth prime be qr, and for a 
logarithm / of length fc, let fr be the restriction of / to the g^-smooth numbers 
between 1 and k. Suppose we have a possibility for fj., with qr ~ exp(^log fc). 
We have therefore used up ^{k,qr) of the elements of as values of / on qr- 
smooth numbers. If we assume suitably random behaviour, the expected number 
of possibilities for fr+i that extend fr is bounded above by 

k{l^^{k,qr)/kY^'''^^^-''^''\ 

since there are ip{k/qr+i, qr+i) numbers between 1 and fc that are ^r+i-smooth 
but not Qr-smooth. Now with u = logfc/logq^ ^ -y/logfc we have tp{k,qr) = 
fcu"""*""^"^ > k^~^ for all large enough fc; and with v — \og{k / qr+i) / y^log k ~ yTogfc 
we have ip{k/qr+i,qr+i) = kv~'"'^°'-^y qr+i > fc^^'^ for all large enough fc. Thus a 
heuristic upper bound for the expected number of possibilities for fr+i (given fr) 
is fc(l — k^'^)'' ' . Since there are 0{k^) = (9(/cCxp(viog fe)^ possibilities for fr, our 
heuristic bound for the number of possibilities for /, summed over all k > ko, tends 
to as fco — > 00. Without some structural reason (such as fc + 1 being prime), this 
suggests that only finitely many fc would achieve Rk as large as exp(-v/log fc) (and 
indeed this bound could heuristically be reduced still further). 

Computations show that Mk = fc for fc < 194, and we computed further that 
Mfe > fc - 1 for fc < 257, but that M258 = 256. 
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